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Abstract. The amount Q of particles that are transported via a path of motion 
is characterized by its expectation value (Q) and by its variance Var(Q). We 
analyze what happens if a particle has two optional paths available to get from 
one site to another site, and in particular what is Var(Q) for the current which 
is induced in a quantum stirring device. It turns out that coherent splitting and 
the stirring effect are intimately related and cannot be understood within the 
framework of the prevailing probabilistic theory. 



1. Introduction 

The possibility to induce DC currents by periodic (AC) modulation of the potential 
is familiar from the context of electronic devices. If an open geometry is concerned, 
it is known as "quantum pumping" [IJ [2l [3l [H [5j [6] , while for closed geometry [3 [8] 
we use the term "quantum stirring" [HI [TU] . Of recent interest is the possibility to 
stir condensed Bose particles that are confined in an optical trap and are described 
by a 3 site Bose-Hubbard Hamiltonian [TTJ [TJJ [T3] . While in a parallel study [T3] we 
explore the role of interactions in this stirring process, in the present work we would 
like to explore a new aspect of the problem that has to do with counting statistics. 
For simplicity we neglect the interactions and accordingly the problem reduces to "one 
particle physics" . 

It is well established [IIlIIiailZlQIlQSllIilE^ that counting statistics 

in the context of shot noise studies provides information on the fluctuations of the 
occupation and on the random probabilistic nature of the quantum-mechanical trans- 
mission/reflection process. In fact these two effects combine together. The prototype 
example is barrier crossing. If the average channel occupation is / and the transmis- 
sion probability is p, then the emerging number of particles Q is proportional to fp, 
while the variance (per-particle) is Var(Q) = (1 — fp)fp. Furthermore the results 
which are derived using classical methods (Master equation; Boltzmann-Langevin) 
are consistent with the quantum mechanical calculation (Scattering approach; Green 
function techniques), and the quantum-classical crossover is related to the statistics 
of the transmission coefficients as reflected by the Fano factor. 

Scope — The purpose of this paper is to argue that the above common wisdom 
does not apply in the context of quantum stirring. In order to demonstrate this point 
we analyze the prototype 3-site system of Fig. 1. We measure the current X that flows 
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through a section in the C\ bond, and define the splitting ratio as A = C\j{c\ + C2). 
The integrated current over a time period is denoted as Q. If the passage probability 
from left-to-right is p we do not get for its variance Var(Q) = (1— Xp)Xp as implied by 
the naive probabilistic considerations, but rather Var(Q) = A 2 (l— p)p. Then we turn 
to analyze a full pumping cycle that consists of two sequential Landau-Zener adiabatic 
passages. During the first half of the cycle A = Aq while during the second half of the 
cycle A = Aq- If Ao 7^ Xq then it follows that the net integrated current is (Q) 7^ 0, 
and we ask what is the variance Var(Q). 

Observations — There are some qualitative observations that are associated with 
our results and we would like to enumerate them in advance: (1) Coherent splitting 
of a wavepackct does not generate a noisy current in the system; (2) The "fractional 
charge" of a partial wavepacket is determined by the splitting ratio, and is physically 
meaningful. (3) The splitting ratio can be greater than unity or negative. This has the 
interpretation of having an induced circulating current in the system. (4) The splitting 
ratio concept allows an easier better understanding of quantum stirring, in comparison 
with the complicated Kubo formalism [7]. (5) The fluctuations in Q reflect the non- 
adiabaticity of the driving cycle. (6) The fluctuations of the integrated current grow 
with time as t and not as \ft. (7) Interference appear differently in the calculation of 
counting statistics when compared with the calculation of occupation statistics. 



2. Definitions 



The current operator I is a conventional observable in quantum mechanics. For a 
single mode ring the current through a section r = ro can be expressed using the 
position and the velocity operators: 

1 = vS(f-r )\ (1) 

We denote by Q the total number of particles that go through the specified section, 
and accordingly define a counting operator: 

Q = fl{t')dt' (2) 







Calligraphic letters are used to distinguish the X and the Q operators, and I(t) = 
U{t)'IUif), where U(t) is the evolution operator. Formally we can diagonalize Q, 
find its Q eigenvalues, and the associated \Q) eigenstates. The initial state of the 
system can be expanded in the Q representation, and accordingly the final probability 
distribution of Q at time t is P(Q) = |(Q|*>| 2 = |*q| 2 - This distribution can be 
characterized by its moments. Of particular interest are the expectation value and 
the variance: 

<Q> = Ei*qI 2 ^ = (*lfi(*)l*> ( 3 ) 

Q 

Var(Q) = (Q 2 )-(Q) 2 (4) 

The naive definition of P(Q) that we have introduced above follows Ref.[25j which was 
later criticized. A more careful analysis l 26 : , 27] of the continuous measurement scheme 
implies that the full counting statistics is characterized by the following physically 
measurable quasi-distribution 

P (Q) = J_ J (^[Te-^ Q ] f [Te + ^ Q ]^e~ lQr dr (5) 
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The naive mathematical definition is obtained if we ignore time ordering: 




(6) 



It is easily verified that Po(<9) has the same first and second moments as P(Q). 
Therefore, for the purpose of variance analysis, we refer below to the latter (naive) 
definition. For more details about full counting statistics in the context of closed 
geometries see [28] . 

3. Modeling 

We consider the simple models that are illustrated in Fig. 1. The Hamiltonian of the 
3 site system is 



Without loss of generality we choose the units of time such as to have the hopping 
amplitude per unit time between |1) and |2) equal 1. For c\ — c-2. — there are 
two energy levels E± = ±1, that correspond to the states |1) ± |2), and Eq = u that 
corresponds to |0). 

The prototype scenario that we consider later is as follows: The particle is initially 
prepared in the left site, and the potential u is slowly raised from u < 1 to u > 1. This 
induces an adiabatic crossing from the left site to the right side. For completeness we 
note that in the standard Bose- Hubbard Hamiltonian all the hopping amplitudes are 
negative, and accordingly the symmetric level E + is in fact the ground state of the 
double well, to which the condensed particles are transported from the Eo = u level. 

More generally we consider later a quantum stirring scenario. The transport is 
measured via the Oh- >1 bond, and accordingly the current operator is represented by 
the matrix 



Assuming control over the couplings c\ and ci we can design a pumping cycle with 
non-zero net transport (Q) ^ 0. See Fig. 2 for illustration. 

4. The two level approximation 

Let us assume a driving cycle of period i p , such that in any moment | ci | , | C2 1 <C 1 and 
u ~ 1, as in Fig. 2. In the strict adiabatic limit the particle stays in the same adiabatic 
level with no leakage to the other levels. But if the driving is not strictly adiabatic 
there is some leakage. Below we introduce the conditions for neglecting the leakage 
to E^. Accordingly the dynamics is very well described within the framework of a 
two level approximation: 




(7) 




(8) 




(9) 



where 



c 




effective coupling 



(10) 
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and 

A = — — — = splitting ratio (11) 

ci + c 2 

If we had only two sites as in the upper illustration of Fig. 1, then A would not emerge. 

One can estimate the transition probability from an initial adiabatic level E n 
to some other adiabatic level E m by writing the Hamiltonian in the adiabatic 
basis: if one changes a parameter X the coupling between the adiabatic levels is 
X [i{dTL/ dX) nm / {E n —E m )\. Then using leading order perturbation theory with 
respect to the driving rate (1/< P oc X), and assuming smooth temporal variation of 
the potentials, one typically obtains (see e.g. Eq.(l3Tj)) 

2 



P 



i „ / 1 
1 



(12) 



where /() provides the time dependence of the coupling between the adiabatic levels, 
and the dynamical phase $(t) is obtained by integrating over E n (X(t))—E m (X(t)). 
The parameter Q characterizes the rate of change of the dynamics phase. In the 
absence of avoided crossings Q is simply the mean level spacing and we label the result 
as Pfgr- This notation implies that we deal with an off-resonant Fermi-golden- rule 
(FGR) transition. But if we have the avoided crossing of the Eq and the E + levels, then 
the predominant contribution to the integral comes from a time interval i LZ = c/ii, 
during which <I> ~ c. Thus for the so-called Landau-Zener transition we effectively 
have in Eq. (|12[) the replacements i — > c and t p i— ► c/u, and we get P LZ P FG r- The 
exact result including the correct prefactors is [29j [30l [31] 

„2" 



P LZ = exp 



u 



(13) 



which is known as the Landau-Zener transition probability. It is important to 
realize that the first order calculation reproduces correctly the singular exponential 
dependence of P on the rate (it) of the driving. 

In this paper we assume an adiabatic stirring process and analyze it within the 
framework of a two level approximation. This means that two inequalities have to be 
satisfied: 

Pfgr « Plz « 1 (14) 

From the above discussion it should be clear that for smooth temporal variation of 
the potential, the first inequality is automatically satisfied simply because c < 1 is 
much smaller compared with the mean level spacing, and it can be further improved 
if we care for separation of time scales t p 3> t L z- 



5. Single path crossing 



Consider a 2 site system as in the upper illustration of Fig. 1. Initially the particle is 
prepared in the left site. Then after some time we measure the probability p to find 
the particle in the right side. In our terminology the probability p characterizes the 
occupation statistics at the end of the dynamical scenario. We assume nothing about 
the dynamical scenario, except of being coherent. This means that the dynamics 
is generated by a 2 x 2 Hamiltonian that can be possibly time dependent. We ask 
whether the counting statistics is related to p. Within a classical probabilistic point of 
view the answer is very simple: one would expect to measure Q = 1 with probability 
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P(l) = p, and Q = with probability P(0) = 1— p. Accordingly one would expect 
to have (Q k ) = p, where k = 0, 1, 2, .... In particular the expectation value would be 
(Q) = P and the variance would be Var(Q) = (1— p)p. 

In the quantum mechanical context this innocent reasoning is wrong. As discussed 



where k = 0, 1, 2, 3, ... and [---J stands for the integer part (i.e. rounded downwards). 
Still this result coincides with the corresponding classical result for k — 1,2, which 
we call restricted quantum-classical correspondence. The purpose of the following 
sections is to explore what happens to the relation between counting statistics and 
occupation statistics once multiple path geometries are involved, and specifically to 
evaluate the first two moments for a quantum stirring process. 

6. Double path crossing 

Consider a 3 site system as in the lower illustration of Fig. 1. Initially the particle 
is prepared in the left site. Then we raise slowly the potential from u < 1 to u > 1. 
At the end of the process there is some probability p to find the particle in the right 
side. Our interest is in the counting statistics of the current that flows during this 
process through one bond, say the i— > 1 bond. Since we have a double path geometry 
we expect a splitting ratio that reflects the relative coupling strength | Ci | 2 / 1 C2 1 2 - 
Using a probabilistic point of view the splitting ratio should combine with transition 
probability (p t— > Xp) , leading to 



Thus for a complete transition (p = 1) with two equally probable paths (A = 1/2) we 
expects to measure the maximal variance Var(Q) = 1/4. 

In the quantum mechanical context this innocent reasoning is wrong. By 
inspection of Eq. © we see that the results for the double path crossing can be obtained 
from the results for the single path crossing using A as scaling factor. Namely, 



The splitting ratio reflects so-to-say the relative weight of C\ in the net hopping 
amplitude ci+C2- For c\ = c 2 we have A = 1/2, which would imply an exact splitting 
of the wavepacket into two equal pieces. In particular we realize that for p — 100% 
transfer efficiency the variance in such a case would not be 1/4 but zero. The value 1/4 
would correspond to a probabilistic (rather than exact) splitting of the wavepacket. 
We may say that the correct description of the transition from the left site to the right 
side is not "the particle has an equal probability to go either via the first or via the 
second path" , but rather "the particle goes simultaneously via the two paths" . Both 
phrasings are equivalent if we have in mind expectation values, but only the latter 
phrasing has the correct connotation once counting statistics is considered. 

Naively we expect that a fraction < A < 1 would be transported via the t—> 1 
bond, while the complementary fraction < (1— A) < 1 would be transported via the 
0n2 bond. But if c\ and c 2 have opposite signs then (say) A becomes larger than 
unity, while (1— A) is negative. This reflects that the driving induces a circulating 




Var(Q) 



(1 - Xp)Xp 



(16) 



(Q) = 

Var(Q) =A 2 (l-p)p 



(17) 
(18) 
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current within the ring, and illuminates the fallacy of the classical peristaltic point of 
view which we discuss below. 

7. The peristaltic picture 

Typically the driving is periodic, and Q is defined as the amount of particles that 
are transported per period. The feasibility to have a non-zero (Q) (non zero "DC" 
current) due to periodic ("AC") driving is known in the context of open geometry 
as "quantum pumping" [3J 0] . We use the term "quantum stirring" [TUl [5] in order 
to describe the analogous effect with regard to a closed device. During an adiabatic 
pumping cycle a conventional two barrier quantum device takes an electron from the 
left lead and ejects it to the right lead. Hence the pumped charge per cycle for a leaky 
pump is naively expected to be (Q) < 1. This naive result is indeed valid if the pump 
is operated in an open geometry between two reservoirs. 

Inspired by the peristaltic picture we assume control over the on-site potential u 
and the coupling constants (ci and C2) which are like "valves". In the first half of the 
cycle C2 = and u is raised across u ~ 1, so as to have an adiabatic passage from the 
left side to the right side via the 1 — >• 1 bond. In the second half of the cycle ci = 
and u is lowered so as to have an adiabatic passage from the right side back to the 
left side via the 1 — > 2 bond. The net effect is to pump one particle per cycle. 

8. Quantum stirring 

We can use the results that have been obtained for a double path adiabatic passage 
in order to illuminate the fallacy of the peristaltic picture once quantum stirring in 
a closed geometry is considered. If during the first half of the cycle A = Xq, and 
during the second half of the cycle A = A5 then in leading approximation (neglecting 
non-adiabatic effects) we get 

(Q) = A -A [per cycle] (19) 

An optional way to derive this result is to make a full 3 level calculation using the 
Kubo formula |7] . Here we have bypassed the "long derivation" by making a reduction 
to an "adiabatic passage" problem. 

It is correct that for a simple minded cycle, where either c\ or ci are zero at each 
stage, we get an agreement with the peristaltic picture. But in general this is a fallacy. 
In fact the essence of quantum stirring is the circulating current which is induced by 
the driving. Contrary to the naive expectation we can get Q > 1 per cycle. This 
happens if C\ and C2 are roughly opposite in sign and hence |A| 3> 1. This also can 
happen if c\ and C2 have the same sign: if we had considered an adiabatic passage at 
u ~ — 1 from |0) to the antisymmetric state — 12) then it would be like replacing C2 
by — C2. In general it is better to say that Q> 1 reflects a huge circulating current 
which is induced if the pumping cycle encircles a degeneracy [7] . 

To avoid confusion it should be emphasized that we are not talking here 
about persistent currents which are zero order conservative effect, but about "linear 
response" which is a first order effect that might have in general both geometric and 
dissipative aspects. It is important to remember that the amount of pumped charge 
per cycle is well defined in the u — > adiabatic limit, implying that it does not depend 
on the actual duration of the cycle. 
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9. Fluctuations 

Having determined (Q) per cycle we would like to find out what is the variance Var(Q). 
The straightforward calculation procedure is to write the current operator in the 
Heisenberg picture: 

I{t) nm = (n\U(t^lU(t)\m) (20) 

and then to integrate it over time so as to get 

-Q\\ iQ± 



-iQl -Q\\ > (21) 

The first two moments (Q) and (Q 2 ) are obtained from this matrix, leading to the 
identifications: 

<Q> = Q\\ (22) 
Var(Q) = \Q ± \ 2 (23) 

For a single path Landau-Zener crossing in a two-site system it has been argued in 
Ref.pS] that the first two moments should be the same as in the classical calculation. 
This is the restricted quantum-classical correspondence that has been mentioned in 
Section 5. Consequently one deduces that 

Q\\ = 1-Plz (24) 



Q± = \/(1--Plz)-Plz x PhaseFactor (25) 

But we have a multiple path geometry, for which restricted quantum-classical 
correspondence cannot be established. Therefore we have to carry out the 
straightforward calculation, which is much more complicated. In practice, in order 
to make the calculation manageable, we can rely on the adiabatic approximation 
scheme. Within this framework the evolution operator is 



U(t) W Y eX P ~* / E n{t')dt' 

n ' L J to 

and accordingly the time dependent current operator is: 
I(t) nm « (n(t)\l\m(t)) x exp 



ra(to) (26) 



/ E nm (t')dt' 

J to 



= M ,„;••,..> " v ) (27. 



If we use the zero order adiabatic states (in u) we get for the diagonal terms zero, 
because the zero order adiabatic states are time-reversal symmetric and hence support 
zero current. If we use the first order adiabatic states we get for the diagonal terms 
a non-zero result with Q» = — Ay. The details of this "Kubo" calculation are not 
included here because this result is a-priory expected on the basis of the much simpler 
analysis of the previous section. 

As emphasized in the previous paragraph, in the case of multiple path geometry 
Qj_ is not related to Qn, and therefore an actual calculation should be carried out. 
The good news is that we get from Eq. (|27|) a non-zero leading order result already in 
the zero order approximation: 

Xce^dt (28) 
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where 



1) 



(2c) 2 alt' 



(29) 



For a single path Landau-Zener transition in a two site system one substitutes A = 1 
and u = id. Then it is possible to demonstrate, see Ref. [28j . that the outcome of the 
integral is ~ %/P L z in agreement with Eq. ([25|) . But we have a multiple path geometry 
for which the result is not known a-priori, so we have to stick to the integral and 
see what comes out. If we had only one Landau Zener crossing we would get \/P LZ 
multiplied by the splitting ratio A. A full pumping cycle is a sequence of two Landau 
Zener crossing, one at t = t\ and the second at t — t 2 . Therefore the integral is a sum 
of two terms, and we get 

2 



Var(Q) 



IQj 



+A e' 



VP2 



P, 



(30) 



where tpi = $(ti) and ip 2 = ®(t 2 ). The result depends on the phase difference 
(p = if2 — fx, which is determined by the time separation of the two crossings. 

For sake of comparison one should realize that the probability p to have remnants 
of the particle in the right side is determined by the coherent addition of the two 
Landau-Zener transitions 1321: 



m/2c 



e^dt 



2 




2 








P LZ 


(31) 



2i-»l + W2c) 2 

Here the interference is with opposite sign, because u in the integrand changes sign. 
Thus we see in a tangible way why due to interference Q± is in general not trivially 
related to Q\\. 



10. Long time limit 

From the above analysis it follows that (Q) grows linearly with the number of cycles. 
An equivalent statement is that the eigenvalues Q± of the Q operator grow linearly 
with the number of cycles. It follows that for a general preparation also the spreading 
•^/Var(Q) grows linearly with time. Needless to say that the probabilistic point of view 
would predict oc \fi growth of the spreading, on the basis of the central limit theorem. 
Thus we have here again another manifestation of the way in which quantum coherent 
behavior differs from its classical stochastic analog. If we have good control over the 
preparation we can select the initial state to be a Floque eigenstate of the quantum 
evolution operator. For such preparation the linear growth of the spreading is avoided, 
and it oscillates around a residual value. 



11. Fractional "charge" 

The derivation of Eq. (fT5|) is based on the observation that the eigenvalues of the 
counting operator are Q = ±^/p. Exactly the same fractional value has surfaced 
in the pioneering publication about counting statistics [25], where the authors had 
interpreted it as an effective "fractional charge" . Their observation was immersed in 
complicated diagrammatic calculations involving a many-body system of Fermions in 
an open geometry. In fact their result has been largely ignored once realized [26] that 
the naive definition of P(Q) is of no physical relevance, while Pq(Q) gives no indication 
for "quantized" fractional charges. 
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It is therefore amusing to realize that a similar idea might emerge in the present 
context. From the above analysis of coherent splitting it follows that for p = 100% 
transfer efficiency it is feasible to a "fractional charge" (Q) = A with (formally) zero 
dispersion. In fact the measured fraction can be greater than unity which we can 
call "mega charge" or it can be negative. The observation of a "mega" charge in this 
context simply reflects the presence of an induced circulating current, which is the 
essence of the quantum stirring effect. 

In any case it should be clear that the notion of "fractional charge" in the context 
of coherent quantum stirring is possibly misleading, and we have raised it merely for 
argumentative purpose. 

12. Summary 

Counting statistics in closed geometries is not of classical nature. Even in the simplest 
problem of a coherent transition between 2 sites, the full counting statistics comes 
out different compared with the probabilistic expectation. Still, the variance comes 
out the same which can be regarded an example of restricted quantum-classical 
correspondence. In contrast to that multiple path geometries require further reasoning, 
because there is no simple way to deduce the counting statistics. 

One observes that the correct description of a quantum passage in a double path 
geometry is not "the particle has an equal probability to go either via the first or 
via the second path" , but rather "the particle goes simultaneously via the two paths" . 
Both phrasings are equivalent if we have in mind expectation values, but only the latter 
has the correct connotation once counting statistics is considered, leading to Eq. fTg]) 
rather than Eq. (|16[) . The coherent splitting of a quantum particle is "exact" rather 
than probabilistic. Furthermore, in a double path adiabatic passage one may find that 
(say) 170% of the particle goes via one path, while —70% goes via the second path. 
This reflects the emergence of a circulating current which is induced by the driving. 

The analysis of adiabatic passage has opened the way to figure out what is the 
counting statistics in the quantum stirring problem. We argue that both the average 
and the spreading of Q grow linearly in time, where the rate of the former characterizes 
the pumping cycle, while the rate of latter depends on the quality of the preparation. 

The result Eq. (f50|) that we have obtained for Var(Q) is exiting because 
it demonstrates how interference gets into the counting statistics calculation 
once multiple path geometries are concerned. Unlike in the calculation of the 
transition probability, the interference is with a different sign, and consequently the 
counting statistics becomes unrelated to the occupation statistics. This should be 
contrasted with the single path crossing problem where non-trivial topology is not 
involved, and consequently the two types of statistics are a-priori related. 

We believe that the analysis of counting statistics in closed geometries that posses 
non-trivial topology not only opens an interesting direction in the study of quantum 
stirring, but also unmasks some essential physics of the counting statistics problem in 
general. 
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Fig.l. Toy models that are analyzed in the Letter: a particle in a 2 site system 
(upper illustration), and a particle in a 3 site system (lower illustration). Initially 
the particle is prepared in the |0) site where it has a potential energy u. The 
hopping amplitudes per unit time (the cs) are also indicated. In the case of a 
3 site system, the time units are chosen such that the hopping amplitude per unit 
time between |1) and \2) equals unity, while the other amplitudes are assumed 
to be small (|ci|, |ca | <C 1). The current is measured through the dotted section. 
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Fig. 2. The adiabatic levels of the 3-site Hamiltonian during one period of a 
pumping cycle. In the absence of coupling (ci = C2 = 0) the E$ = u(t) level 
intersects the symmetric E+ = 1 level. With non zero coupling these intersections 
become avoided crossings, and the particle follows adiabatically the thickened 
lines. For presentation purpose we indicate that either ci or C2 equal zero 
("blocked"), but in the general analysis we allow any splitting ratio, including 
the possibility ci = C2 of having the same amplitude to take either of the two 
paths. 



